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Abstract 

Recent analytical and numerical work on field driven domain wall propagation in nanowires has 
shown that for large transverse anisotropy and sufficiently large applied fields the Walker profile 
becomes unstable before the breakdown field, giving way to a slower stationary domain wall. We 
perform an asymptotic expansion of the Landau Lifshitz Gilbert equation for large transverse 
magnetic anisotropy and show that the asymptotic dynamics reproduces this behavior. At low 
applied field the speed increases linearly with the field and the profile is the classic Landau profile. 
Beyond a critical value of the applied field the domain wall slows down. The appearance of a slower 
domain wall profile in the asymptotic dynamics is due to a transition from a pushed to a pulled 
front of a reaction diffusion equation. 
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Magnetic domain wall propagation is an active area of research both as an interesting phys¬ 
ical phenomenon as well as for its possible applications in logic devices, magnetic memory 


elements and others [l|. The dynamics of magnetic domain walls is described by the Lan- 
dan Lifshitz Gilbert (LLG) eqnation j2, 3| which cannot be solved analytically except in 
very special cases. For an infinite medinm with nniaxial anisotropy and an external field 
applied along the symmetry axis, the Walker solntion jd] provides the best known analytical 
expression for the profile and speed of the domain wall. The exact Walker solntion, predicts 
that the speed increases linearly with the field np to a critical field Hw Above this valne 
a sndden drop in velocity and an irregular precessing motion of the magnetization appears. 
Field induced domain wall propagation in thin films and nanowires has been examined with 

n 

greater detail in recent work. The numerical study showed that depending on the relative 
magnitude of hard axis anisotropy different scenarios arise. For small hard axis anisotropy 
the Walker solution is realized. For sufficiently large values of the hard axis anisotropy the 
Walker breakdown does not occur. There is a slowdown of the domain wall due to spin 
wave emission and no sudden drop in speed. For the largest values considered in the 
domain wall speed changes from a regime of linear growth with the applied field to a regime 
of slower growth with increasing applied field. This last behavior is observed when both 
the exchange constant and uniaxial anisotropy are much smaller than hard axis anisotropy. 


Further numerical studies 


, ITI] analyze in detail the nature of the spin waves emitted and 


distinguish two scenarios, depending on the relative values of the exchange and anisotropy 


constants. The parameter ranges studied in 


and 


differ, however, in both cases the 


Walker breakdown is not observed when the transverse anisotropy is sufficiently large. The 
stability of the Walker solution with respect to small perturbations has been studied recently 
81] using dynamical systems techniques. The analysis of the spectrum of a perturbation to 
the Walker solution shows that it may become absolutely or convectively unstable before 
the breakdown field. This instability is found for sufficiently large transverse anisotropy and 
for fields larger than a critical value. 

The purpose of this work is to study the dynamics of the LLG equation for a nanowire when 
the transverse anisotropy is large by means of an asymptotic expansion. The asymptotic 
expansion captures the slower relaxation dynamics of the domain wall and filters out the 
fast spin waves js]. We find that the leading order asymptotic dynamics predicts a transition 
from a Walker type regime to a regime with slower domain wall motion. In leading order 
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the dynamics of the in plane magnetization obeys a reaction diffusion equation, and the 
perpendicular magnetization is slaved to the in-plane components. The slowdown of the 
domain wall in this asymptotic regime appears as a transition from a pushed to a pulled 
front at a critical value of the applied held. 

The starting point of the calculation is the LLG equation for the magnetization. The 
material has magnetization M = Msfh where Mg is the saturation magnetization and m = 
(mi, m 2 , m 3 ) is a unit vector along the direction of magnetization. The dynamic evolution 
of the magnetization is governed by the LLG equation. 


dM ^ ^ M dM 

— = -joM X ft. + a— X ^ 


( 1 ) 


where is the effective magnetic held, 70 = \'y\fio, 7 is the gyromagnetic ratio of the 
electron and /iq is the magnetic permeability of vacuum. The constant a > 0 is the dimen¬ 
sionless phenomenological Gilbert damping coefficient. We consider a thin and narrow him 
in the {x,y) plane, with the easy axis along its length. The strip is subject to an applied 
magnetic held along the easy axis Ha = HaX. The him is thin and narrow so that the magne¬ 


tization may be assumed 


lOj to depend on the easy axis coordinate, M{x, y, z) = M[x). In 


addition, in this geometry the demagnetizing held has a local representation as an ehective 


perpendicular anisotropy so that, as in 

Gp. d‘^M 2K 


8|, the ehective magnetic held is given by 


H,g = HaX -f- 


+ 


;MiX 


2K, 


7 M 3 Z, 


( 2 ) 


yoM^ dx‘^ yoMg yoMg 

where Cex is the exchange constant, the easy axis uniaxial anisotropy and the per¬ 
pendicular anisotropy. 

Introducing Mg as unit of magnetic held, and introducing the dimensionless space and time 
variables ^ = X\JKu/Cex and r = yo\'y\Mgt we rewrite equations (II]) and (|2|) in dimensionless 
form 


dm 


dm 


—— = —m X heft + am x —— 


with 


dr 


Kg = haX -F -K\ 


dr 




+ K\\mix — K±m^z. 


(3) 


(4) 


2''" de 

where ha is the dimensionless applied held and the dimensionless numbers that have appeared 
are K\\ = 2Ku/{yoMg), iLj, = 2Kd/{yoMf). Equations (I3|) and (01) describe the dynamics 
of the problem. 
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We are interested in the case of a perpendicular anisotropy much larger than the uniaxial 
in plane anisotropy. In this situation the perpendicular magnetization m 3 will be smaller 
than the in plane components. We will also assume that the dimensionless applied field is 
weak. We search then for a solution of the LLG equation in the asymptotic limit ha « K±, 
<< K± and therefore, mi, m 2 » m 3 . Let then 

mi = mio + emu + • • •, m 2 = m2o + em2i + ..., m 3 = emso + e^m 3 i + ..., 


where e is a small quantity. Since the perpendicular anisotropy is larger than the uniaxial 
anisotropy and the applied held is weak, we introduce the scaling iLy = eK\\, ha = eha with 
K± of order one. The components of the effective magnetic held, hefr = (hi, ^ 2 , h^) become 
then hi = ehm + e^hn + ... with the leading order components given by 


1 ~ (9^m 

hio = ha+ 2'^ll Q^2 —^ ^\\mioi 


ho = and hso =-i^^rngo. (5) 


Furthermore we introduce a slow time scale s = er and notice that the leading order com¬ 
ponents of the in plane magnetization satisfy 

m\o + mlQ = 1 - 0 ( 6 ^). 


Introducing these scalings in Eq.([3]) and expanding in e one obtains 

dmiii 


ds 

dm2o 

ds 


— ~m2oh3o + 0(e), 
= miohgo -|- 0(e), 


0 = -mioh 2 o + m 2 ohio + a{mio—^ - m 2 o-^) + 0(e). 
Substituting fITal) and fl7b|l into fITcD and using ([H]) we hnd that in leading order, 

hgo = — (mioh2o — m2ohio) 


a 


and equations (l7a]l and (!7b|l become 

9m 10 


__ m2o 

ds a 

dm2o mio 


(mioh2o - m2ohio) 


Because of 


„ (mioh2o - m2ohio). 

ds a 

we can write mio = cos 9, m 2 o = sin 9. Using ([5]) in equations 

do 1 ~ 

a— = -K\\9cc — sin 9(ha + K\\ cos 9) 
ds ^ 


mgo = 


1 d9 


K\ ds 


( 6 ) 

(7a) 

(7b) 

(7c) 

( 8 ) 

(9a) 
(9b) 
we obtain 
( 10 ) 
( 11 ) 
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Finally going back to the unsealed time variable r and parameters we write the 

leading order magnetization components as 


mi = cos 9, 7712 = sin 0, m^ = 


do K\\ 

a— = -^9fe — sin 9{ha + K\\ cos 9). 

OT 2 M ^ 


1 dO 
K I dr 


where 


(12a) 

(12b) 


Equations fll2p show that the leading order dynamics is determined by the equation for 
the in-plane magnetization components, the perpendicular magnetization is slaved to the 
tangential magnetization. Equation fll2bp is the well studied reaction diffusion equation, for 


which we know that an ini 
front of minimal speed 


ial perturbation to an unstable state evolves into the monotonic 


12j. In order to render fll2bp into the standard form we introduce 


the dependent variable u dehned hj 9 = 7r(l — u) which satishes 


aUr = + f{u), 


with f{u) = —sm.nu{ha — Ku costtm). 

TT 


(13) 


The diffusion constant D = K\\/2 and the reaction term / satisihes f{u) > 0 in (0,1), 
/(O) = /(I) = 0. A small perturbation to the unstable state m = 0 (0 = tt) evolves into a 


traveling monotonic front of minimal speed c* 


Ini, El 


that joins the unstable state to the 


stable state u = 1{9 = Ob The minimal speed can be obtained from a variational principle 


13l | and is bounded by 


m 


ckpp = -\/Df'{Q) <c* < -\/D sup/(m)/m. 
a a 


(14) 


When the upper and lower bounds coincide the speed is exactly ckpp and the traveling front 
is called a KPP or pulled front. 

In the present problem Eq. flT^ has the exact traveling front solution 

2 

“ arctan 

71 

This solution is not a KPP front, it is a so called pushed front. This is the front into which 
an initial condition will evolve it is effectively the front of minimal speed. It is not difficult 
to verify that as ha increases this is not the speed of the front. For ha > 4iF|| the upper 
and lower bounds in ffTTl) coincide and the speed of the front must be the KPP value. The 
transition from a pushed to a pulled front may occur before the upper and lower bounds 




, ha D 

where cn = —\ 

a \ All 


(16) 


5 

















coincide. In this problem for which there is an exact solution we know that the transition 
will occur when cn = ckpp- That is, 


c = 


hg 

a 


I 


if ha < 2Kii 

:^D{ha-K\\) iiha>2K\\ 

Going back to the physical variables, we have then that the speed of the domain wall is 
given by 

a\f^^ho\h\Ha 


if Ha < 


V = 


- 2Ka 


iiHa> 


IiqMs 


( 16 ) 


olMs V s a u a fiQMs 

In the small held regime Ha < H^ = the magnetization prohle is obtained from 

firsjl and it is given by 


nil = tanh 


(X 


vt) 


A 


m 2 = sech 


(X 


vt) 


A 


iloMsHa 

m 3 = —;=-sech 

V2aKa 


(X 


vt) 


A 


. (17) 


where the domain wall width is given by A = (1/2) For an applied held larger than 

He we cannot construct an explicit solution for the magnetization, we can only determine the 
speed. The general theory of reaction dihusion equations guarantees that it is a monotonic 
decaying front similar in shape to fflSj) . Equations flT^ and flTTl) constitute our main result. 
These results, obtained from the LLG equation in the case << K\ and for a weak applied 
held, explain qualitatively the results of the numerical simulations p-lTj and of the stability 
results | 8 |. At low helds the speed of the front is proportional to the applied held Ha and 
inversely proportional to the damping coefficient a. The magnetization prohle and the speed 
share the main features of the Walker solution, the velocity shows linear dependence on the 
applied held and inverse proportionality on the damping constant a. The Walker breakdown 
held Hw = olK^I2 is of order one, and therefore large compared to the transition held 
He- Thus, we recover the behavior described in 5|-^: for sufficiently large perpendicular 
anisotropy the Walker solution loses stability before the breakdown held to a slower moving 
domain wall. When the applied held is weak and Kn « K± the numerical integrations 


m 


BE 


7| show that the speed increases slowly with the held once the Walker solution loses 


stability, in agreement with the results found in this work. 

The asymptotic approach that we have used is based on , where the numerical simulations, 
(although for a diherent demagnetizing held), show that the asymptotic dynamics reproduces 
the relaxation dynamics of the full LLG equation, hltering out the spin waves. Reaction 
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diffusion dynamics has also been encountered in thin nanotubes 


14| . where the Walker 


breakdown is not observed. A transition from a fast to a slower domain wall also occurs in 


thin nanotubes as reported in 


15 


, ll6|. The asymptotic dynamics of the LLG equation has 


been studied by several authors in different limiting parameter ranges, and wave type motion 


3 , 0 . 


In the present problem we 


governed by other evolution equations has been derived 
have chosen a parameter regime for which recent numerical work has been performed and 
found qualitative agreement with the results reported in them. 
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